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Abstract 

We show that the rise in F% at small x and large Q 2 seen at HERA is indeed the non-Regge 
double asymptotic scaling behaviour expected from the perturbative emission of strongly ordered 
hard gluons. An alternative explanation, in which there is no strong ordering, and a new hard 
Reggeon is generated, is also tried but found wanting: its theoretical short-comings are betrayed 
by its failure to properly account for the HERA data. 



Talk given at "The Heart of the Matter" , Blois, June 1994, 
to be published in the proceedings (Editions Frontieres). 



CERN-TH.7421/94 
September 1994 



On leave from a Royal Society University Research Fellowship. 
On leave from INFN, Sezione di Torino, Italy. 



The structure function F^x^Q 2 ), recently measured |1J by the HERA experiments 
ZEUS and HI in the hitherto unexplored region 10 -4 ^ x ^ 10 -2 , 5GeV 2 & Q 2 & 
10 5 GeV 2 , rises dramatically both as x decreases and Q 2 increases. The form of this rise 
may be clearly exhibited by using the variables 0] 

^\A*¥K' p^ln^/ln^, (1) 

where t = ln(Q 2 /A 2 ). Rescaling F^ by a simple multiplicative factor R' F = Na 1 / 2 pe Sa '< p ', 
where 5 = ||, we may then plot it on a logarithmic scale against a. The resulting plot 
(fig. 1) is interesting in two respects: firstly when a is large enough all the data lie on 
a single line, quite independently of the value of p (provided p too is large enough), and 
secondly the rise of log R'pF^ with a is linear. The slope of the rise0 is 2.37 ± 0.16 
(dotted line in fig. 1). 

Remarkable though it is, such a rise was not unanticipated: the behaviour of F2 
at small x was computed in perturbative QCD soon after the discovery of asymptotic 
freedom[Q, using the operator product expansion (at leading twist), the renormalization 
group, perturbation theory (at one loop) and assuming that at some low starting scale Qo 
the small x behaviour of F2 is given by conventional Regge behaviour, and thus flat. The 
resulting perturbatively generated non-Regge behaviour takes the form 

F 2 ~AT/(J)(^)- 1 e 2 ^, (2) 

where Af and / ~ 1 + 0(p~ 1 ) depend on the (soft) non-perturbative input at Qo- Since 
27 = 4^N c /Po = 2.4 is simply a numerical constant, perturbative QCD thus predicts 
precisely the double scaling behaviour described above: when F 2 is rescaled by a factor 
Rp = e~ 2l ° R' F to remove both the linear rise and the sub-asymptotic effects R' F , the 
resulting structure function is asymptotically independent of both a and p. 

The rise in F2 is thus generated perturbatively through iteration of the simple pro- 
cesses g — > gg and g — > qq: the splitting functions P gg (x) and P qg (x) are both singular 
at leading order as x — > 0. In fact at small x the gluon evolution equation reduces to a 
wave equation in light-cone variables £ = log — , £ = log -p, with a (negative) mass term 
which generates an exponential growth in xg inside the light-cone (see fig. 2: turn it anti- 
clockwise through 45°). This then drives a similar growth in F%. This explains why the 
natural variables to use when x < xq and t > to are a = v^C an d p = VwC) an d also why 
a generic exponential rise (a-scaling) will be produced isotropically (p-scaling) whenever 
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the boundary conditions set at x = xq and t = to are sufficiently soft (and in particular 
g(x;to) ~ x~ x as x — > 0), so that the main source of the waves is close to the origin. 

Double asymptotic scaling seems to set in just as precociously as Bjorken scaling: 
almost all the HERA data falls within the scaling region. For this reason we take Qq = 
1 GeV; higher values would imply large double scaling violations not seen in the data. To 
illustrate this point we display in fig. 3 the MRS prediction |J DO', a soft distribution fitted 
to pre-HERA data, and then evolved (at two loops) from Qq = 4 GeV 2 . Although similar 
in shape to the double scaling curves (shown dotted in fig. 3), DO' agrees less well with 
the HERA data. Dropping the starting scale to 1 GeV would have given a prediction which, 
while having the same shape as double scaling, also had the correct normalization. [|5| 

Note that although 1 GeV is at the boundary of the nonperturbative region, we only 
compare to data with Q 2 ^ 5GeV 2 , where higher twist corrections to F2 are negligible. 
Higher twist evolution below this scale may be absorbed into similar uncertainties in the 
starting distribution. 

While double scaling captures the essential physics of the rise in F2, there are several 
corrections still to be included: 

i) Sub-asymptotic corrections, both at one and two loops; these may (and have been0||) 
computed explicitly; the dominant uncertainty is now the shape of the initial gluon distri- 
bution. 

ii) Post-asymptotic corrections, due to the leading singularities in the anomalous dimen- 
sions at higher orders in a s ; these are now known explicitly and may be included in the 
Altarelli-Parisi equation to all orders. In this way the the strong kx ordering implicit at 
leading order is gradually relaxed in a controlled way as we go to smaller x. In the HERA 
region such corrections are still quite small however H.* Subleading singularities might 
also be calculable. 

iii) Higher twist (screening) corrections, necessary to restore unitarity post-asymptotically; 
these are probably very small in the HERA region 0, though they are difficult to compute 
reliably. 

Thus in the 'soft pomeron' approach (perturbative evolution of an initially soft non- 
perturbative input), higher order corrections are all (except possibly iii) under control, in 
the sense that they are small in the HERA region, and may furthermore be systematically 

* Shift the axes of fig. 9 of ref. [|| by jq and \ respectively, to take into account our choice of 
xq and Qq: almost all the HERA data then lies within these authors' convergence criterion. 
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computed within a well denned scheme. This is essentially because the operator product 
expansion and renormalization group are used to control the expansions. By contrast in 
the 'hard pomeron' approach (see ref. [|IU| for a critical review) these fundamentals are 
abandoned in favour of a particular diagrammatic resummation. This inevitably leads to 
difficulties. 

The starting point for this latter approach is the 'BFKL equation' which resums logs 
of 1/x at fixed Q 2 ; kr ordering is abandoned at the outset. Solution of the equation gives 
asymptotically (at small x and large Q 2 ) the Regge-like behaviour 



F 2 ~MJQ 2 /Q 2 x- X , A = 121og2^. (3) 

v 7T 



Such a rise violates double scaling || , and indeed is much stronger than that seen at HERA 
(it would imply an exponential rather than linear rise in fig. 1, for example). However it 
need not be taken too literally, because a number of important theoretical problems remain 
unresolved: 

I) The BFKL equation does not consider evolution in Q 2 , and thus does not resum leading 
logs. The equation is thus inconsistent with the renormalization group, and its solution 



cannot be related to known distributions at larger x. Alternative equations |TT| [T^], which 
also have no kr ordering but include Q 2 evolution, are difficult to solve, but seem to give 
numerical results consistent with leading order Altarelli-Parisi evolution. This is because 



when the coupling runs, the hard pomeron cut dissolves into a sequence of poles[|T3[, and 
the right-most pole has only a small residue. 

II) When the kr ordering implicit in the conventional leading twist evolution equations 
is abandoned, contributions with different twist become mixed together. This leads to 
nonperturbative infrared effects which cannot be easily factored out, except by projecting 
onto the leading twist. An attempt to solve this problem may be found in ref.[|l4]]. The 
implementation of the (kinematically necessary) ultraviolet cutoff on fcy is also problem- 
atic. UT5J Sub leading corrections, while probably important (to ensure momentum conserva- 
tion, for example), are as yet unknown. 

III) Since unitarity is now violated more quickly since the growth (4) is much stronger 
than (2), screening corrections are now much larger JIBJ ; they are also even more difficult 
to calculate. 

While most of these problems remain unresolved, it is not yet really possible to com- 
pare the 'hard' approach with the data. Despite (or perhaps because of) this, it has 
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recently become popular to employ a 'hybrid' approach, in which the perturbatively based 
small- a; behaviour (3) is used as a nonperturbative boundary condition at Qq, which is then 
evolved using the conventional two loop evolution equations. The resulting Q 2 dependence 
is then relatively weak, as the Lipatov cut is to the right of the leading singularities in 
the anomalous dimensions. fllOfl Since such an approach makes little sense theoretically, it is 



rather satisfying that the data themselves contradict it: a hard starting distribution propa- 
gates inside the light-cone anisotropically, generating strong violations of p-scaling.0This 
is demonstrated rather clearly by the MRS prediction^ D- ; , displayed in fig. 4 (the 'hard' 
pomeron curves of ref.0, also with A = |, are shown dotted for comparison). In fact D-' 
fits the data even less well than DO', essentially because the sharp rise at small x is not 
matched by a corresponding rise at large Q 2 . Of course by reducing A and/or using a 
boundary condition with only a very small relative admixture of hard to soft behaviour, 
it is still possible to fit the data: the recent MRS H fit mimics double scaling in precisely 
this way. * 

In conclusion, the 'hard' pomeron is neither a proven consequence of perturbative 
QCD, nor is it as yet visible in the structure function data. Indeed, the simple prediction 
of perturbative QCD[Q has been beautifully confirmed |2[] @ by the measured[Q rise in F2 
at HERA. This is fortunate, because it means that perturbative QCD can be used to make 
reliable predictions at small x. It should thus be possible to use the HERA data to further 
test QCD, measure ct s , and predict QCD backgrounds at the LHC by extrapolating to yet 
higher values of a. 

Acknowledgements: We would like to thank S. Catani, R.K. Ellis, F. Hautmann, 
Z. Kunszt, J. Kwiecinski, P.V. Landshoff, E.M. Levin, R.G. Roberts and D.A. Ross for 
discussions. 



* In fact there is only around 3% of the hard solution in MRS(H): see eqn(24) of ref. [17] 
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Figure Captions 

Fig. 1. The rise in F%: R'pF^ vs - a - Also shown is the double scaling prediction||, and 
a fitted straight line for comparison. 

Fig. 2. The x-t plane, showing the different kinematic regions, and in particular the 
double scaling region, with scaling co-ordinates a and p. 

Fig. 3. Double scaling plots of RfF^ vs. a) a and b) p. The data are taken from ref. 0, 
and the curves are those of the MRS parton distributions DO'. The double scaling 
curves [0 are shown dotted. 

Fig. 4. As fig. 3, but with the curves now corresponding to D-' and the 'hard pomeron' 
curves computed in ref. (dotted) . 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9409373vl 



This figure "fig2-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9409373vl 




Q ~1GeV 



This figure "fig2-2.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9409373vl 



This figure "fig2-3.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9409373vl 



This figure "fig2-4.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9409373vl 
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